Geotropic tracers in turbulent flows: a proxy for fluid acceleration 
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Abstract 

We investigate the statistics of orientation of small, neutrally buoyant, spherical tracers whose 
center of mass is displaced from the geometrical center. If appropriate-sized particles are considered, 
a linear relation can be derived between the horizontal components of the orientation vector and 
the same components of acceleration. Direct numerical simulations are carried out, showing that 
such relation can be used to reconstruct the statistics of acceleration fluctuations up to the order 
of the gravitational acceleration. Based on such results, we suggest a novel method for the local 
experimental measurement of accelerations in turbulent flows. 



1 



I. INTRODUCTION 



The Lagrangian investigation of turbulence has dramatically improved in the last few 
years in experimental techniques, theoretical models and numerical simulations 1]. These 
progresses benefited from the increased range of Reynolds numbers accessible for investiga- 
tion (in particular, for simulations) and the improved accuracy of measurement techniques. 
On the theoretical side, we have now phenomenological models able to quantitatively explain 
;he Lagrangian properties of turbulence such as, e.g. the statistics of velocity increments 
^ and accelerations Grounded on the successes of Lagrangian investigations, recent 
experimental and numerical studies started to investigate the motion of complex objects 
in turbulent flows (J-^. The motivations are both fundamental and applicative. In this 
short note we suggest a possible technique to measure turbulent accelerations without the 
need of particle tracking, by means of the local measurement of the orientation of finite-size 
particles. The idea relies on spherical particles whose average density is that of the carrier 
fluid (so that they are neutrally buoyant), but whose center of mass is displaced with respect 
to the geometrical center (implying that the orientation is determined by the gravitational 
torque and that due to the fluid). By means of direct numerical simulations (DNS) we show 
that information on particle orientation can be used to estimate fluid accelerations up to 
the order of gravitational acceleration. 

The paper is organized as follows. In Sect. [Ill we discuss the theoretical basis of the 
technique in its simplest implementation. Sect. II III presents some preliminary validation of 
the method based on numerical simulations. Finally, Sect. [IV] is devoted to conclusions and 
perspectives. 



II. THE MOTION OF GEOTROPIC TRACERS 



We consider the trajectory, x(t), of a neutrally-buoyant sphere small enough such that 
its dynamics can be approximated by that of a passive tracer, 

-^ = u(x,t), (1) 

transported by a flow u(x,t). As sketched in Fig. [TJ we assume that the particle center 
of mass C is displaced by a distance h with respect to the geometrical center O (which 
is the center of buoyancy). The displacement determines the particle orientation, defined 
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FIG. 1: An example of geotropic particle with black and white pattern. 



by the versor p directed opposite to the center of mass. In a steady fluid, the equilibrium 
orientation points towards the vertical direction z opposite to gravity. In the presence of 
a flow, the direction p is determined by the balance between the different torques acting 
on the particle. Because of particle asymmetry, an external force /, such as gravity mg, 
exerts a torque Tf — hp x f. In addition, the particle immersed in a fluid experiences a 
viscous torque T v = 87ri>pr 3 (u}/2 — ft), where a> = V x u is the fluid vorticity, v and p 
are the fluid kinematic viscosity and density, respectively, and Q is the angular velocity of 
the sphere. If the particle Reynolds number is very small, we can assume creeping flow 
conditions around the sphere, which impose equilibrium between the external forces and the 
viscous ones, in this particular case zero total torque Tf + T V = 0. Considering that we can 
write ft — (Q ■ jiyp + p x p, from p x (Tf + T v ) = we end with the following equation for 



the orientation 



9] 



^ = -—[A 

dt 2v 1 



[A-p)p] + -ujxp, 



(2) 



where A denotes the acceleration on the sphere and the constant Vq = 3u/h, having the 
dimension of a velocity, weights the contribution of external forces to particle orientation. 
We remark that in the case of axisymmetric non-spherical particles an additional term is 
present in ([2]) [l 9]. 

In a fluid at rest, the only external force entering Eq. ([2]) is gravity, so that A = g = 



(0, 0, —g). The result is that particles orient upwards in a time 0(vq/ g). Such phenomenon 
is well known in bio-fluid-dynamics, and is at the basis of the ability of some bottom-heavy 
ohytoplankters to swim towards the sea surface (a phenomenon dubbed negative gravitaxis 
qI), maximizing the exposition to light and thus the photosynthetic activity. 

In a turbulent flow, advected particles are subject to intense accelerations, so that, locally, 
gravity must be corrected due to inertial forces. The total acceleration A acting on the 
particle is thus given by A = g — a, where a = du/dt is the fluid acceleration at the particle 
position. 

In the following, for the sake of simplicity, we assume that fluid acceleration is smaller 
than gravity, i.e. g 3> a rms . As we will see this is true for flows at moderate Reynolds 
numbers only, but this assumption greatly simplifies the analysis of 02]). Formally, we can 
write A = g — ea in (J2J), with e a small non-dimensional number. Further we consider 
the limit of fast reorientation i.e. v — > (notice that this limit can be controlled by 
the displacement h as, vq = 3u/h) and we neglect the viscous term u> x p in (j2j), which is 
consistent when vou rms <C a rms . With these approximations, Eq. (j2J) reduces to A = (A-p)p, 
which explicitly reads 

ea x = (ea-p + gp z )p x 

ea y = (ea-p + gp z )p y (3) 
ea z + g = (ea ■ p + gp z )p z . 

From (jHJ) one can see that the orientation vector must have the form p = (eq x ,eq y ,l), 
indeed since p 2 = 1 the correction to p z will be 0(e 2 ) so that we can neglect it at this 
level. Plugging the expression for p in Eq. (jHJ), at 0(e) we obtain q = (a x /g,a y /g,0). In 
conclusion, measuring the orientation of the particle with respect to the vertical we can 
measure two components of the fluid acceleration 

This result is valid under the assumption that a rms g and, therefore, in general for not 
too high Re. This condition together with the smallness of the vorticity term implies fast 
orientation (vo/g t v ). As a consequence, if vq is too large, the first effect we expect is 
that vorticity becomes relevant, with an increase of the tilting angle, so that using fll]) could 
lead to an over-estimate of accelerations. In more general conditions, it is in principle still 



possible to use (|2J) to gather information on the acceleration statistics but this requires less 
direct procedures, which we will not consider in this preliminary study. 



III. NUMERICAL SIMULATIONS OF GEOTROPIC TRACERS IN TURBU- 
LENCE 

In this section we illustrate the behavior of geotropic particles in realistic turbulent flows 
and explore the range of validity of the result fll]) by means of DNS of the dynamics of 
geotropic tracers together with the Navier-Stokes equations for an incompressible flow. Tra- 
jectories (up to 2xl0 5 ) are stored together with a and uj in statistically stationary conditions. 
Equation (T5]) is then integrated starting from random orientations and for different values of 
Vq. After an initial transient of the order of Vo/g, during which particles forget their initial 
orientation, we can compare the acceleration a with the prediction of Eq. 

We have performed simulations of homogeneous-isotropic turbulence by means of a paral- 
lel pseudo-spectral code in a cubic box with periodic boundary conditions at Re A ~ 200 with 
resolution 512 3 . Statistical stationarity was maintained via a Gaussian, delta-correlated in 
time, random forcing at small wave-numbers. Eq. (JT|) is integrated evaluating the velocity 
at particle position by means of trilinear interpolation. Moreover, we have also exploited 
a database [lOj of previously simulated Lagrangian trajectories at resolution 2048 3 and 
ReA — 400, for which acceleration and vorticity was available. Equation (J2J) was integrated 
using a second-order Adams-Bashforth scheme. In order to get physical relevance from the 
DNS we rescale space and time with dimensional values. This is easily done by matching 
Kolmogorov scale and time with experimental results at similar ReA, as shown in Table [B 
We remark that this rescaling is not unique as ReA fixes a ratio of scales (and times) and 
not an absolute scale. This point is crucial as the parameter vq is limited by the size of the 
particle and g is obviously fixed. In the following we will use laboratory experiments with an 



integral scale of the order of few cm for rescaling our simulations to physical values 



Fig. [2] shows an example of time series of the two components of the acceleration a x and 
a y obtained following a Lagrangian tracer in the flow at ReA = 200. The initial condition of 
the orientation is along the z axis, p(0) = (0, 0, 1). The red line represents the acceleration 
obtained according to (pO) from the x component of the orientation vector, a x = gp x of a 
particle with vq ~ 0.006ms -1 . The corresponding relaxation time under gravity is r = 
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Re A 


rj (m) 




u' (ms 1 ) 


e (m 2 s 3 ) 


L (m) 


a' (ms 2 ) 


r(a) 


200 


191 x 10~ 6 


37 x 10~ 3 


0.037 


7.1 x 10~ 4 


0.07 


0.14 


8.22 


400 


98 x 10~ 6 


9.5 x 10~ 3 


0.097 


9.1 x 10~ 3 


0.10 


0.93 


2.58 



TABLE I: Parameters of the simulations made dimensional on the basis of laboratory experiments 



at similar Reynolds numbers 
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11]. rj = (i^ 3 /e) 1//4 is the Kolmogorov scale, t v = (vje) 1 ! 2 the 
Kolmogorov time, v! is the root mean square of the velocity, e is the energy dissipation for unit 
mass.L = u' 3 /e is the integral length scale, a' = e 3 / 4 zy -1 / 4 is the root mean square acceleration. T 
is the integration time. For both simulations g = 9.8ms -2 . 

Vo/g ~6x lCT 4 s. In this case a rms <C g and, therefore, the estimation (jlj is fully justified 
and indeed the acceleration is reproduced quite accurately. 

In Fig. 13] we show an example for a trajectory in a turbulent flow at Re\ = 400. Although 
the rms of acceleration a rms ~ g 3 / 4 // -1 / 4 is smaller than g, particle experience fluctuations 
comparable to, or even larger than g, where the assumptions leading to (jlj) are not applicable. 
These large fluctuations of Lagrangian acceleration are typical in turbulence and physically 



correspond to event of trapping of tracers in small scale vortices [3|, ]12|, ]13j. As shown in 
Fig. [HJ during these events the orientation vector p is unable to fully follow the acceleration 
fluctuation, which results to be slightly underestimated. 

On a more quantitative level, Fig. H] shows the probability density function (PDF) of 
acceleration compared with the estimation obtained via (j3J). For each value of ReA considered 
we simulate the results of three hypothetical experiments, with particles of different sizes. 
As discussed above, geotropic orientation is expected to be a good proxy for acceleration 
only in the limit of small Vq, i.e. for fast orientation. As apparent from both panels in 
Fig. HJ when a large enough displacement h is considered, the statistics of acceleration is 
reproduced remarkably well by particle orientation. However, this is not the case if less 
biased particles are considered. This clearly implies a lower limit in the size of particles 
used, a factor that must be taken into account in the design of possible experiments. As 
mentioned above, vorticity can be neglected only if VoU rms /a rms ~ vq/5u v < 1. By applying 
the definition of the Kolmogorov scale du^rj/u = 1 and that of vq, the constraint reduces (a 
part from order-one coefficients) to rj ^ h. This inequality can pose a problem both for the 
validity of §B) and the actual statistics seen by the particle. Both points will be discussed in 
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the final section. As for now we will just consider this condition in the frame of our model, 
assuming that the corrections are small as long as the particle size is of the same order as r\. 

If (as in our case) one considers a set of experiments all using water and with comparable 
integral scales, an increase in Re corresponds to a smaller viscous scale, thus decreasing 
the minimum particle size required to reconstruct acceleration. As an example of this we 
considered the case of particles with h = 0.2mm. As evident from FigJH using Eq. (JJJ) on 
statistics obtained with such particles would lead to an overestimate of larger accelerations 
at ReA = 200 (triangles in FiglJJi) while they would be acceptable candidates at ReA = 400 
(triangles in panel b). However, the largest acceleration that can be measured by means of (J2J) 
is g. For experiments at higher Re where very large accelerations are present, this introduces 
a cut-off in the estimated accelerations. As evident from the results at ReA = 400 the core of 
the PDF is approximately correct, even if values above 0.5 -r-0.7g are underrepresented. We 
stress that the simulations exhibited accelerations up to 80 a rms (not shown for graphical 
reasons), while g ~ 6.3a rms if rescaled over the experimental parameters. Analysis of the 
variance of acceleration performed for ReA = 200 (inset of Fig. Hfci) reveals that the second 
moment of the distribution is correctly recovered asymptotically in h/rj. The same cannot 
be verified at ReA = 400, since the cut-off at g prevents convergence of the second moment 
of estimated accelerations. 

IV. CONCLUSION AND DISCUSSION 

Summarizing our numerical results, it appears that the orientation of biased particles 
could be a viable proxy for fluid acceleration, at least at moderate Re. It is clearly important 
to establish a way to estimate the proper particle size based on the parameters of the 
turbulent flow to be examined. 

Although particle size does not enter directly in the model, the offset h clearly puts a 
lower limit on particle radius. This point must be carefully considered. Particles should be 
sufficiently biased to ensure dominance of acceleration over viscous rotation due to vorticity, 
but too large particles would not obey the assumptions leading to (jSJ). 

On the other hand, experimental limitations should be considered. In order to use 01]) to 
directly measure fluid acceleration one has to measure the tilt angle of a geotropic particle 
transported by the flow. One possibility is to use small spherical particles with the upper 
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and lower hemisphere of different color as in the example of Fig. [TJ By measuring the angle 
9 of the particle "equator" with respect the horizontal plane one has p x = sin8. 

A precise determination of 9 requires sufficient resolution of the particle pattern and 
therefore not too small particles. On the other hand, because the measure is instantaneous, 
there is no need to follow the particles. Therefore, the camera could be placed to zoom a 
small region of the fluid only, and to acquire data when a particle comes in that region. 

Let us finally comment about possible corrections to the described behavior, for the two 
cases of particles too small or too large. If the offset h is too small, the vorticity term in 
(|2J) is no longer negligible. As a consequence, reconstruction of acceleration would require 
independent information on vorticity, so that a more complex method would be required. 
Furthermore, fast orientation is at the basis of (jl]) which allows one to avoid particle tracking, 
and would be important to follow high frequency fluctuations accurately. In the case of too 
large particles, the creeping flow assumption would be inaccurate. Nonetheless, the tilting 
angle would still provide information on the fluid acceleration but equation (j2j) has to be 
modified to take inertial terms into account. 

A further aspect that should be taken into account is finite size effects on particle trajec- 
tory. In general one expects that particles larger than the Kolmogorov scale deviate from 
fluid trajectories. However there is evidence that acceleration statistics are not strongly 



influenced by particle size 



14j-|l6|. Numerical and experimental results suggest that addi- 



tion of the so called Faxen terms in the equation for particle trajectory can account for the 
main deviations, providing a method to estimate the related errors JjJIla]. Such corrections 
should become relevant when the radius of the particle is larger than r)\/Re\, which for 
experiments comparable to the ones we considered would give O(10)^ 16| thus allowing for 
some range of sizes to explore. 

In spite of the above discussed limitations, we think that the idea of exploiting biased 
particles to measure acceleration without tracking may be interesting especially if technology 
can be pushed to the possibility to measure the tilting angle of many particles at the same 
time, allowing for the reconstruction of the spatial field of accelerations. 
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FIG. 2: x component (upper plot) and y component (lower plot) of the acceleration of one particle 
computed from numerical simulations at Re^ = 200 (black line) together with the acceleration 
estimated from the x and y component of the orientation vector, i.e. a x = gp x a y = gp y see dU, of 
a geotropic particle with vq = 6mms _1 corresponding to a displacement h = 0.5mm. 
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FIG. 4: PDFs of acceleration in one horizontal direction for ReA = 200 (a) and ReA = 400 (b). 
Estimates obtained according to §2$ (X = gp x , symbols) are compared with fluid acceleration 
{X = a x , line). Three values of particle bias were used, which rescaled on experimental values 
correspond to h = 0.1mm (circles), 0.2mm (triangles) and 0.5mm (squares). By comparison, it 
is evident that the intermediate value gives a good estimate at higher ReA but is not satisfying 
at the lower one(see text). In panel (b) the of g (vertical lines) marks the upper cutoff for 

measurable accelerations. . In the inset of panel (a): relative error on the estimate of a = J (of) 



